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ABSTRACT 

I have undertaken a comprehensive statistical investigation of the ultracool 
dwarf companion distribution (spectral type M6 and later). Utilizing a Bayesian 
algorithm, I tested models of the companion distribution against data from an 
extensive set of space and ground-based imaging observations of nearby ultracool 
dwarfs. My main conclusions are fivefold: 1) Confirm that the concentration of 
high mass ratio ultracool binary systems is a fundamental feature of the com- 
panion distribution, not an observational or selection bias; 2) Determine that the 
wide (> ~20 AU) binary frequency can be no more the 1 — 2%; 3) Show that the 
decreasing binary frequency with later spectral types is a real trend; 4) Demon- 
strate that a large population of currently undetected low mass ratio systems are 
not consistent with the current data; 5) Find that the population of spectroscopic 
binaries must be be at least 30% that of currently known ultracool binaries. The 
best fit value for the overall M6 and later binary frequency is ~20% — 22%, of 
which only ~6% consists of currently undetected companions with separations 
less than 1 AU. If this is correct, then the upper limit of the ultracool binary 
population discovered to date is ~75%. I find that the numerical simulation re- 
sults of the ejection formation method are inconsistent with the outcome of this 
analysis. However, dynamics do seem to play an important role as simulations of 
small-N clusters and triple system decays produce results similar to those of this 
work. The observational efforts required to improve these constraints are shown 
to be primarily large spectroscopic binary surveys and improved high-resolution 
imaging techniques. 



Subject headings: stars: low- mass, brown dwarfs; (stars:) binaries: visual; stars: 
formation; methods: statistical 
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Introduction 



An understanding of the star fo rmation process has long been a major goal of modern 
astronomy. The pioneering work of ISafpeterl (119551 ) established the Initial Mass Function 
(IMF), the number of stars born per unit mass per unit volume, as a fundamental distri- 
bution for star forma tion theory. The study of the IMF was later expanded by the work of 
Miller fc Scald (119791 ) and countless efforts since in many young clusters and the field. A 
number of functional forms of the IMF have been tested in varied environments, such as the 



lognormal (jMiller fc Scalolll979l ) and the power-law (jSalpeterlll955l ). As a community, we are 
convinced that there is a link between the IMF and the global process by which those stars 
formed. However, there is another set of by-products of star formation that can elucidate 
the formation process: binary companions. 

The frequency and distribution of phy sical properties of companio ns to intermediate to 
low-mass stars has been well characterized. iDuquennoy fc Mayorl (Il99ll ) sh owed that G stars 
have a binary rate of ~60% with a Gaussian distribution of separations. iFischer &: Marcy 
( 119921 ) performed similar studies but for early-M dwarfs and found a binary rate of ~40% 
with a similar Gaussian distribution in separation to the G dwarf studies. Both of these 
studies used samples drawn from catalogs of known objects that have their own selection 
biases. They also used multiple observational techniques to derive their binary frequencies 
and companion distributions. As a result, the frequencies derived incorporate corrections 
to account for the inherent uncertainties in the original samples and their combination. 
This is an unfortunate, but necessary step in most companion analyzes (including this one). 
Hotter stars binary frequencies are not as well understood. This is mostly because there 
are relatively few very massive stars and they evolve very quickly. However, there is a well- 
established trend that the binary frequency increase with mass at the high end. For example, 
the b inary frequency of O stars is at least 75% and most likely much higher (IMason et al 
1998h . 



Recent years have ushered in a new area of binary studies: very low- mass stars (< 
O.15M ), which includes ultracool dwarfs (spectral type M7 and later), many of the coolest 
and lowest mass extending into the brown dwarf regime. After years of fruitless searching 
for brown d warfs , the profusion of deep far-red optical (Sloan Digital Sky Survey, SPSS 



York et al.l (120001 )) and near- infrared (Deep Near I nfrared Survey, DENIS lEpchtein et al. 



(119971 ). and the Two-Micron All Sky Survey, 2MASS ISkrutskie et al.l (120061 )) sky surveys in 
the last decade catalyzed the discovery of hundreds of nearby ultracool dwarfs. This enables 
us to probe an entirely new mass range of the star formation process. 

Shortly after the discovery of this large population of near by ultracool dwa r fs in the 
field, the community began to search for companions to them, see iBurgasser et al.l (120071 ) or 
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Duchene et al.l (120071 ) for a review. The observed field ultracool dwarf companion distribution 
has been shown to be quite different from that of higher mass stars. These imaging surveys 
have shown that ~15% of late-M and L dwarfs are resolved as binary systems. Only a 
handful of these systems has a separation exceeding 15 AU. This is in contrast to ultracool 
dwarfs with higher mass primaries; for example, VB 10, the archetypal l ate-type M dwarf, 
lies 400 AU from its primary, the M3 dwar f, Gl 752A (Ivan Biesbroeck!ll944f). while the nearby 



T dwarfs, Gl 229B (INakaiima et al.lll995l ). Gl 570D (IBureasser et al.ll2000h . and e Indi Bab 
(IScholz et al.l 120031 ) . are all wide components in multiple systems. On the other hand, a 
handful of wide (> 50 AU) , predominantly young (<10 Myr), ultracool binaries have been 
identified. iLuhmanl (l2004bl ) has identified a pair of late-type M dwarfs separate d by 240 AU 



i n p O p hiuchus, an additiona l 212 AU M dwarf pair in Ophiuchus was found by lClose et al. 



(120071 ). IChauvin et al.l (120041 ) have discovered a very low- mass brown dwarf companio n of 
the TW Hya member 2M1207— 39, with a separation of 60 AU, and ILuhmanl (120061 ) has 
found a 210 AU pairing in Chameleon I. There are two field binary systems discovered so 



far: D ENIS 0551-4434AB a M8.5/L0 pair with a separation of over 200 AU (IBilleres et al. 
2005 ): and 2 MASS 0126AB, a M6.5/M8 pai r with a separation of 5100 AU (Urtigau et al. 
2003). It should be noted that lArtigau et all (120071 ) argue that 2MASS 0126AB may be part 
of the TW Hydra association, and thus may be a young, wide binary like 2MASS 1207-39. 

The detected field binary systems are predominantly equal brightness/mass systems 
(Figure [1]) in tightly bound orbits (Figure [2]) with relatively few exceptions. However, we 
would like to know if these features only due to observational/selection biases or are real 
features of the ultracool dwarf companion distribution. For example, the majority of the 
imaging surveys from which these conclusions are drawn are based on magnitude limited 
samples. Such samples tend to pick out the brightest objects of a given spectral class. Equal 
brightness/mass binaries are a prime contaminant because they appear to be very bright 
when the are in fact two separate objects. In any analysis of the overall distribution of 
ultracool dwarf binaries such biases must be taken into account. In this work, I will use 



techniques based on those derived in lBurgasser et al.l (120031 ). 



Numerical simulations of the various formation scenarios are just beginning to pro- 
duce prediction s of th e binary properties of samples formed via those mechanisms, see 



Burgasser et al.l (120071 ) for a review of formation mechanisms and simulation results. These 



simulations encompass a wide variety, including standard cloud fragmentation, massive cir- 



outcome of my analysis to the few simulations that produce testable results 


(Bate, Bonnell, & Bromm 


2002, 


2003: 


Bate & Bonnell 


2005; 


Sterzik & Durisen 


2003; 


Umbreit et al. 


2005 


)• 



To date there has not yet been a comprehensive statistical analysis of the ultracool 
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dwarf companion distribution. This work will provide the first such study for companions to 
nearby field late-M, L, and T dwarfs, using only imaging data. I use a Bayesian algorithm to 
determine the underlying companion distribution, utilizing both the detections in the surveys 
as well as the non-detections. The Bayesian method also allows the easy incorporation of 
multiple data sets with one another to generate a unified conclusion. 

In this paper I will bring together the disparate ultracool bina ry imaging data sets 



with theoretical models that describe the evolution of brown dwarfs (IBurrows et al.ll200ll ) 



to constrain t he ultracool binary distribution via a Bayesian statistical approach similar to 



that found in lAUen et al.l (120051 ). hereafter referred to as A05. I will also determine if the 
unusual features in the ultracool binary distribution (the paucity of wide systems and the 
preference for near-equal masses) are real or are an observational/selection bias. This paper 
is organized as follows: Section 2 discusses the binary surveys to be used in my statistical 
analysis. Section 3 explains my basic Bayesian methodology and Section 4 outlines how 
I produce the binary distribution models. Section 5 discusses the results and Section 6 
summarizes my conclusions. 



2. Summary of Companion Surveys and Data Sets 
2.1. Surveys and Detections 



As will be seen in Section HJ the Bayesian algorithm requires data sets that report 
a complete list of discovered binaries as well as non-detections. This eliminates systems 
discovered serendipitously from this analysis. Additionally, I will focus this work on present 
imaging results, thus excluding spectroscopic binaries. The reasons for this are twofold: 1) 
The spectroscopic surveys are not yet as fully developed as the imaging surveys in size and 
scope; 2) A great deal of additional modeling is needed to properly include the spectroscopic 
results and it was deemed outside the reach of this paper. The analysis of spectroscopic 
results will be t he topic of a future pap er. Two examples of ultracool binaries to be excluded 
are PP1 15AB (IBasri &: Martini Il999l ) and elndiBab, the former was discovered via high 
resolution spectral measurements of a brown dwarf candidate in the Pleiades and the latter 
serendipitously found in a proper motion survey. The serendipitous binaries cannot be 
included because the Bayesian algorithm I use requires a full set of detections and non- 
detections. Such data cannot be obtained from a chance discovery. In all ~20% of known 
ultracool dwarf binary systems fall into t he excluded category. This is based on the compiled 
list of 75 ultracool binaries in Table 1 of iBurgasser et al.l (120071 ) and eliminating those found 
either serendipitously or in non-imaging surveys. 
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The majority of the ultracool dwarf comp anion survey s have been carried out with HST. 
T he first set, u sing WFPC2, is summari zed in lBouv et al. (|2003h . and is the combined work 
of iReid et al.l (120011 ) , iGizis et al.l (120031 ) , and iBouy et al.l (120031 ) . These gr oups surveyed a 



total of 131 late-M and L dwarfs with HST and found 25 binary systems. iBurgasser et al. 



( 120031 ) also used WFPC2 to examine 10 T dwarfs and resolved 2 binaries. Additional HST 
surveys, using the NICMOS NIC1 camera, obtained hi gh resolution ima ges of 3 M dwarfs and 
49 L dwarfs tha t are m ostly within 20 pc of the Sun (IReid et al.l 120061 ) . and 22 T dwarfs by 
Burgasser et al. ( 2006h. These observations resolved a further 14 binaries. The other high- 



resolution surveys (IClose et al.ll2003l ; ISiegler et al.ll2003l . 120051 ) are a ground-b ased program 



that u sed the nea r -IR a daptive optics system at Gemini- North, Hokupa'a (IGraves et al. 



19981). IClose et al 



(200 3J) exa mined 39 M8 through L0 dwarfs and found 9 binary systems, 



and ISiegler et al.l (120031 . 120051 ) examined 41 M6-M7.5 dwarfs and resolved 5 systems. 



This analysis also uses data from two low-resolution, wide-field, ground-based surveys 
of late-M and L dwarfs for com panions. The preliminary results of a su rvey using the Keck 
Teles cope (IKoerner et al. Il999h . with the Near InfraRed Camera (NIRC) (IMatthews fc Soiferl 
19941 ). found 3 binary systems out of ten surveyed. The full data set (76 objects) is not 
yet published (Koerne r et al., in prep.) and includes one other binary system, which was 
previously resolved by IReid et al.l (120011 ). All four of these binary systems have also been 
observed within the HST programs described above. The final survey is one that used 
NSFCam on NASA's IRTF t o study 133 late-M and L dwarfs and found no candidate 
companions (jAUen et al.l 120071 ). 



In all, a total 361 individual targets were imaged with these surveys, including 53 that 
are resolved as binary systems, which yields an observed binary fraction of ~15%. Figure 
[3] displays the spectral type and distance distributions of the primary star population for 
the above surveys. These systems represent a large cross-section of the late-M., L, and T 
dwarf populations at the time of each survey, all with s election biases. The early surveys 



Koerner et all Il999l: llleid et all I2OO1I ; IGizis et al.l I2OO3I : iBouv et al.l I2OO3I : IBurgasser eiTa! 



2003 



2006 



Allen et alil2007h observed all the object s they knew about a t the time reg a rdless 



of dis tance or selection effects. The AO surveys of IClose et al.l (120031 ); ISiegler et al.l (12003 



2005|) were constrained to the brightest and closest objects based on the Hokupa'a natural 



gu ide star AO syste m they were using. The most recently devised and implemented survey 
of IReid et al.l (120061 ) has a much better defined sample of ultracool dwarfs. These objects 
were primarily taken from those found in lCruz et al.l (120031 ). and represents a volume-limited 
sample of early and mid L dwarfs within 20 pc of the Sun. 



The binaries found thus far point to a companion distribution that differs from that 
of higher mass primaries, G, K, and early-M stars. Figure [1] displays a comparison of the 
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mass ratio distribution of the detected binaries in several nearby companion populations. 
The trend in these data, continued by the ultracool dwarf companions, is to an increasingly 
peaked mass ratio distribution with decreasing primary mass. It remains possible, though 
unlikely , that this is du e the o bservational bias towards equal brightness binaries. Recent 
work of iBurgasser et al.l (120061 ) demonstrates that, when this bias is taken into account, the 
observed ultracool binary frequency is considerably lower than that of higher mass primaries 
when compared in similar regions of parameter space. Thus, this trend appears to be genuine. 



2.2. Window Functions and Bias Correction 



The surveys included here are not sensitive to low-mass ratio companions at small sep- 
arations and they suffer from decreasing linear resolution with increasing primary distance. 
This can be seen in Figure HI which displays the observed window functions of all the surveys 
used in this analysis. I define the window function as the area in observational space (delta 
magnitude vs. the log of the projected linear separation) in which a companion could have 
been detected. This was determined for each individual target based on information in their 
respective source papers. For the HST WFPC2 sur veys of late-M and L dwarfs I used the 
typical sensitivity limits given in iGizis et al.l (120031). and the WFPC2 surveys of T dwarfs 
used the individ ual sensitivity est imates given in IBurgasser et al.l (120031 ). The HST NICMOS 
L dwarf survey (IReid et al.ll2006l ) provided a sensitivity curve as a function of radial distance 
from the cent ral source and was us ed to prepare that window function. The NICMOS T 
dwarf survey (IBurgasser et al.l 120070 again provided individual sensitivity m easurements for 



each target. For the AO surveys (I Close et a 



20031 : ISierier et al.112003 . l2005h . I used the sen 



sitivity curve provided in lSiegler et al.l (120051 ) for all the observations, as the author indicated 
that the curve was representative (Nick Siegler, private communication). The Keck survey 
also contains individual sensitivity measurements for each target and they were used in the 
construction of t heir window funct ions (David Koerner, private communication). Finally, 



the IRTF survey (jAUen et al.l 120071 ) also provides individual sensitivity measurements. 



The window functions displayed in Figure H] are constructed as follows. In each bin of 
observable space (delta mag vs. projected separation) I examine the sensitivity curves de- 
scribed above for the appropriate survey to which a target belongs. If that bin is observable, 
it is assigned a 1, if not it is assigned a 0. The only exception to this assignment is near the 
edges of each field. All of these surveys were conducted with square fields of view (FOV). 
However, the projected separation is not probed uniformly out to the value represented at 
the corner of the field. Consequently, I created a weighting function that assumes a circular 
FOV out to a radius equal to half the diagonal of the square FOV which scales the area 
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probed outside of a radius equal to the FOV/2. This yields a gradual drop off of sensi- 
tivity to larger projected separations. The final modification of the window function is to 
correct for the selection bias of the samples from which the targets were compiled. Recall, 
that magnitude limited surveys are biased towards the brightest objects, which are more 
likely to be unr e solved binary systems. To correct for this, I use a technique described in 
Burgasser et al.l (120031 ). They derive the following form to correct this bias as a function of 
delta magnitude: 

(1 + io°- 4xAm ) L5 (1) 

The result of this correction is to decrease the significance of detections at small delta 
magnitudes, or near-equal mass systems, as required. 

The general structure in the window functions is dominated by a step-like shape at small 
separations. This is the result of decreasing sensitivity to low mass ratio binaries the closer 
you get to the central source. The cutoff at a delta magnitude of 5 in the HST WFPC2 
window function is due to the average sensitivity curve used for each observation in that 
survey. This is different from the K-band window function which is a composite of four dif- 



feren t surveys flClose et al.ll2003l : ISiegler et al J 12003 . l2005t lAUen et all 120071 : Koerner et al. 



19991 ) . These surveys were able to probe to smaller mass ratios (larger differences between 
primary and secondary mass) because smaller companions are brighter in the near-IR pho- 
tometry ba n ds th an in the optical WFPC2 bands. Also, both Koerner et al. (in prep) and 
Allen et al.l (120071 ) provide individual sensitivity measurements for each target which leads 
to a non-uniform cutoff at large delta magnitudes. Overall, there are two major holes in 
the coverage of these window functions: 1) very wide, low mass ratio companions and 2) 
sub-AU separation companions. As a result, there may be a population of tight or wide, low- 
mass companions hidden where they cannot be seen in the current surveys. These potential 
populations are discussed in Section [51 



3. Models of the Substellar Companion Distribution 

The first step in the Bayesian statistical analysis is to build a physical companion 
distribution model that I will test against the data. There are two main components to 
consider in a companion distribution, how they are distributed around each star (orbital 
parameters) and their relative masses (mass ratio). Here I setup and justify those model 
distributions. 
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3.1. Separation Distribution 



The separation distribution of ultracool dwarf binaries will be modeled as a Gaussian. 
This is based on existing data for companions to GKM stars, as described in Sections 1 
and 2. The projected separation distribution of companions to GKM stars were found to 
be consistent with a log- normal (see Figure 7 of IDuquennoy &: Mayorl (ll99ll ) , Figure lb of 
Mayor et al.l (119921 ). and Figure 2b of IFischer fc Marcvl Jl992h) . Figu re [2] displays the sepa- 
rations of ultraco ol binaries, as summar i zed in iBurgasser et al.l (120071 ). superimposed on the 
G star binaries of iDuquennoy fc Mayorl (119911 ). While the sma ll separations are incom plete, 
the overall observed distribution has been shown in the work of IBurgasser et al.l (120071 ) to be 
fit well by a symmetric Gaussian. This is similar to the separation distrib utions derived for 
higher mass primaries (IDuquennoy fc Mayorl Il99ll ; IFischer fc Marcvl Il992f). However , I will 
also consider an asymmetric Gaussian shape, as described in I Jeffries fc Maxtedl (120051 ). They 
use this shape to statistically examine the likelihood of spectroscopic binaries at separations 
less than 1-2 AU and introduce a sharp cutoff at 10-20 AU to account for the lack of wide 
companions. In interest of completeness, I will test this model as well. 

The parameterized model form of the standard symmetric Gaussian is straightforward: 

1 



P(log(o)| log(a ), a a )Alog(a) 



(log(a)-log(q ))" 
2^ 



2ira n 



Alog(a) 



(2) 



where log(a) is the log of the semi-major axis (AU) of the orbit, Alog(a) is the bin size in 
separation, log(a ) defines the center point of the log-normal sep a ration model, and a a is 
the width. The asymmetric Gaussian model of I Jeffries fc Maxtedl (120051 ) is constructed by 
multiplying Equation [2] by a strong exponential cutoff beginning at 10 AU. 



P(\og(a) | log(ao), (T a ) Alog(a) 



2no. 



(log(a)-log(a )) _ 

. e 2^ e (i-iog(a)) 3 Alog(a) 



(3) 



However, the semi-major axis is not what is measured in a binary system; it is the 
projected separation on the plane of the sky. In order to transform the above model, the 
companions are randomly distributed in inclination, phase, and eccentricity. The projection 
onto the plane of the sky for a particular value of the inclination angle, i, the phase angle, 
<ft, and the eccentricity, e, is given by the following form: 



cos 2 (0)sin 2 (i) + sin 2 (0) 



, ,. ,., v . (4) 

1 + ecos(0) V VYV v ' VYV v ' 

where s is the projected separation in AU. When distributed over all values of i (0° — > 90°), 
<p (0 —>■ 27r), and e (0 — > 1) this transforms a single value of a to a range in s. Figure [5] 
displays a comparison of the original and the transformed separation distribution models for 
both the symmetric and asymmetric Gaussian models. 
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3.2. Mass Ratio Distribution 

Recall Figure [U which displays the known mass ratio distributions for G, K, early-M 
stars, and the currently known ultracool binary systems. Also recall that the mass ratios 
of G and K binary systems are widely distributed, with a peak near q = 0.3. This begins 
to change as one looks at lower mass primaries. Binaries amo ng the early-M dw arfs and 



the local solar neighborhood, which is dominated by M dwarfs ( jReid fc Gizid 119971 ). have a 
mostly-flat mass ratio distribution Therefore, there seems to be a general evolution toward 
high mass ratio binary systems with decreasing primary mass. The ultracool dwarf mass 
ratio distribution appears to continue this trend. 

A power law is used to describe the mass distribution of brown dwarf binary systems. 
However, for the binary distribution it is not the mass that is interesting, but the mass ratio, 
q = MijM\. So, the mass ratio distribution is assumed to follow a power law, g 7 , where 7 
is the power law index. This power-law is assumed to be constant and does not allow for 
the roll over or cutoff of the binary mass ratio distribution. Thus it may overestimate the 
number of small mass ratio companions.. The index 7 is not exactly similar to the field or 
young cluster IMF power-law index a, because 7 is for the combined binary system mass 
ratio not an individual mass. Thus, if both the primary and secondary are drawn from a 
similar IMF then the slope of the mass ratio function (7) will be close to 0. As stated earlier, 
this is true for intermediate mass G stars. However, if 7 is found to be close to zero than it 
can be said that the low mass binary distribution is drawn from a similar mass function as 
that of single low-mass stars, (~ — 0.3, A05). 

The overall normalization of the models is factored in at this point. The normalization, 
N, is set to be the overall binary fraction, the percentage of ultracool dwarfs that have lower 
mass companions. When combined with the above power-law in q and a minimum q value 
of 0.02, the following form is obtained: 

1.0 

P{q\N, 1 )Aq = q^AqxN/ ( 5 ) 

g=0.02 

where Aq is the bin size in mass ratio, and N/ Y2l=o 02 ^ 7 ^ s ^ ne normalization. The left 
panel of Figure [6] displays model mass ratio distributions for two values of 7, —1 and +1. 
The 7 = — 1 model predicts relatively few equal mass systems and an increasing number of 
unequal mass systems. However, for 7 = +1 the distribution is opposite that of 7 = — 1, 
with more equal mass systems than unequal. 

These models then need to be transformed to the observable, the magnitude difference 
of the secondary from the primary at a particular band pass. The techniques use d to trans- 
form these models are nearly identical to those used in I Allen et al.l (120031 . 120051 ). I use the 
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Burrows et al.1 (120011 ) ultracool dwarf evolutionary models to provide a link between mass, 
bolometric luminosity, effective temperature, and age. However, what is measured is not 
the bolometric luminosity or e ffective surface temperat ure but a magnitude in a particular 



bandpass. As in A05, I use the iGolimowski et al.1 (120041 ) bolometric corrections to transform 



the Burrows models into the observable magnitudes from the bolometri c. The bolometric 



correc tions are assigned based on the effective temperatures derived in IGolimowski et al. 



( 120041 ) . The main difference between A05 and this work is that I am now concerned with 
a mass ratio distribution and not plain mass. To determine these ratios the transformation 
requires information on the primary stars. As stated in A05, each late-M, L, and T spectral 
type has a unique distribution in mass and age. This is important because the spectral type 
of each potential primary is known. With the range of spectral types sampled in a particular 
survey, a distribution of mass and age for each primary type and the corresponding range of 
mass and age for possible secondaries can be constructed. 

Since M6 is the earliest spectral type considered in this analysis, a maximum mass of 
0.15 M & is considered. The minimum mass is a bit trickier to determine. Substellar objects 
have a mass-age degeneracy. Put simply, as a brown dwarf of a particular mass ages it also 
cools, thus its spectral type changes. The latest spectral type in this analysis is a T8. To take 
into account this floor to the analysis, I take the absolute magnitude of T8s in the different 
band passes. This maximum magnitude translates into a range of masses depending on the 
age. As in A05, I assume a flat age distribution. Thus the mass-age degeneracy of brown 
dwarfs is accounted for and is built into the model companion distributions. 

The right-hand panel of Figure [6] displays the mass ratio distributions transformed to 
Arrix for 7 equal to —1 and +1. For 7 = — 1, the structure in the resultant distribution 
is similar to the luminosity functions generated in A05 Figures 2 & 3. This is because, 
with a power law index of —1, the distribution of companions will follow a similar form as 
single field dwarfs do because that index yields an increasing number of objects with small 
q. However, the 7 = +1 distribution is strikingly different. This model predicts very few 
unequal mass systems (large A magnitudes), and results in a solitary peak at Arrix = 0, 
or equal mass systems. These models are combined with the separation distribution models 
from the previous subsection to form a complete set of companion distribution models for 
the Bayesian analysis of the companion distribution. 



4. Bayesian Analysis 



Before I can proceed with my statistical analysis of the ultracool binary surveys of 
Section [2] to determine the underlying companion distribution, we need an understanding of 
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my statistical methodology. This section provides a brief tutorial on the Bayesian method 
and how I apply it to constrain the ultracool binary distribution models. The methods 
described here are similar to those in A05. 



4.1. Bayesian Core Algorithm 

The core of this method is Bayes' rule: 

Pm\{O})ccP({O}\{0})xP({9}) (6) 

where 6 is the model, {0} is the set of observed data, P({#}|{0}), the posterior distribution, 
is the probability of the model given the data, P({0}|{#}), the likelihood function, is the 
likelihood of the data gi yen the mo del, and P({9}), the prior distribution, is the initial 



probability of the model (lSivialll996l ). The output posterior distribution provides a wealth 
of information on the model parameters, from the best fit parameter set to correlations 
between parameters. Using Bayes' rule I calculate the probability that these data would 
have been measured given a hypothesized model (the likelihood function). However, I wish 
to know the probability that a hypothesized model is true given the measured data (the 
posterior distribution). The power of Bayes' rule lies in the simple relation of these two 
quantities. 

The specification of a prior distribution remains the most controversial aspect of Bayesian 
analyzes and must be considered carefully. The prior folds previous observational and the- 
oretical evidence into the analysis in more than one manner. Prior information can also be 
included in the functional forms taken by the models, e.g. the Gaussian shape of the sepa- 
ration distribution. One technique uses the posterior distribution from a previous analysis. 
This enables the same analysis to be performed in light of new and improved data. In this 
way, one can iterate over multiple data sets, thereby incorporating them into the analysis of 
a single set of models to provide one unified result. This is ideal for the study of the ultracool 
dwarf companion distribution because there is no single data set for all such dwarfs, and new 
data are continually made available. Priors can also be constructed if no previous knowledge 
of the problem exists. In this case the prior distribution should not impart a bias on any 
parameter value under consideration. These types of priors fall under the broad heading of 
conjugate distributions. Since there has been relatively little work done on the statistical 
analysis of the ultracool dwarf companion distribution, the initial prior distributions used 
here will be constructed using conjugate distributions. What little prior information there 
is went into the selection of the forms of the model distributions. 
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4.2. Basic Problem and Setup 

The Bayesian analysis of the companion distribution is carried out in the same way as it 
was for the field mass function in A05. A Poisson form is chosen as the likelihood function. 
Recall the general form of the Poisson likelihood function: 

In this case, the model observations, as derived in Section [31 include R, the number of 
companions per star as a function of Am and log(s). R is defined as R({0}, {0}), where 
{0} represents the observable space (Am, log(s)), and {8} is the vector of model parameters 
(N, 7, log(ao), (To) as defined previously in Section [3j 

N det and N obs must also be defined and constructed. They are simply two dimensional 
distributions of the number of detections, N det ({0}), and the number of times a bin in ob- 
servable space is sampled, N o f, s ({0}). N o f, s ({0}) also describes the window function, in other 
words, a measure of the areas of observational space covered and the number of observations 
each area has with respect to the others, see Section 2.2. With all three components of the 
Poisson likelihood function defined, it takes on the following form at one point {{0}ij) in 
observable space: 



P({O} hJ \{0}) oc R({O} t ,\{0}) Nd ^ {O} ^xe- 



■Nobsaoh^xRaohjio}) 



where the subscript i represents the ith bin is log(s), projected separation, and j the jth bin 
in Am, delta magnitude. 

The priors on the model companion distributions must also be determined. The shape 
the prior takes depends on the functional form of the models. What previous knowledge of 
the ultracool companion distribution there is has gone into the derivation of the functional 
forms chosen in Section [3j I wish to fully test the entire range of parameters permitted 
by these models, thus, I will choose the conjugate priors associated with the model forms 
I have assumed. Th e conjugate distributions for the parameters of a Gaussian distribution 



are well-established ( lSivialll996l ). In this case, the prior on the lognormal center separation 
parameter, log(ao), is flat or constant. The width, a a , is a scale parameter. A maximum 
entropy argument can be used to determine the prior distribution on a a . To do so, the 
constraints to which the distribution on a a is subject must be defined. The only constraint 
on the probability distribution of a a is an invariance to changes in scale; i.e., the units 
of distance can be changed from AU to centimeters with no effect on the outcome. This 
means that a a is a scale p arameter, and the most 'ignorant' prior distribution is given by 



P(cr a ) oc l/<7 a (jSivial Il996h . The normalization, N, is also a scale parameter; it is required 



to be independent of the units of the calculation. The power-law index, 7, of the mass ratio 
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distribution is unconstrained and does not depend on the scale, so it is given a flat prior. 
The complete prior distribution is given as follows: 

^(W)oc^xl (9) 
N a a 

Equations [HJ and [9] together form the RHS of the unnormalized Bayes' rule (Equation 
[6]). This form is only true for one point in observable space, {0}ij. However, it can easily be 
manipulated to be both easier to calculate and cover the entire observable space. Assuming 
that the observation of the companion rate in an observed bin, {0}i j, is independent of 
the measurement in any other bin, the likelihoods and priors at each observed point can be 
multiplied to obtain the following form: 

p({6}\{0}) « p(W)x n^({0}| W) M (10) 

where the subscripts are the same as in Equation [HJ The prior distribution is outside the 
product because it is a function of the model parameters and not the observable dimensions. 
The calculation becomes easier, in the same way as the field mass function analysis in A05, 
if the natural logarithm is taken of both sides of Equation [TUJ It takes on the following form: 

ln(P({6}\{0})) = ln{P{{6})) +J^N deti .({0})xln(R({0}\{9}) i , j ) - 

#ote w ({O})xfl({0}|{0}) W + constant (11) 

where the sum is over the two observable dimensions, Am; and log(s)j. This is the final 
form of the likelihood function, and it will be used in Section [5] to constrain the binary 
distributions models. 



5. Results 

To perform the Bayesian calculations, the models, described in Section [31 and the data, 
from Section [21 are fed into Equation [TH The output posterior distribution is four dimen- 
sional (log(a ), a a , N, 7). As a result, it is impossible to display the entire distribution, so 
it is marginalized over different sets of parameters. This process simply means the collapse 
of the posterior distribution along different parameter axes. The marginalized distributions 
of each individual parameter are displayed in Figures [7] and [HJ, for the symmetric and asym- 
metric Gaussian models respectively. The best fit parameter values, with la uncertainties, 
are shown in Table 1. Each individual parameter distribution is well-behaved and exhibits 
a clear peak. 
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5.1. Model Parameter Results 

The parameter values returned by the two different models are very similar (Table 
1). This is not surprising because the region where they differ (separations >10 AU) has 
very few detections. However, we should recall that many of the detections that have been 
made at these wide separations have been excluded from this analysis. Both parameter 
sets are consistent in finding a peak to the separation distribution near 7 AU, an overall 
binary frequency in the low 20%'s, and a steep mass ratio spectrum (7 ~ 1.8). The main 
difference lies in the Gaussian width parameter, a a . The best fit values are quite similar 
but the la uncertainties are different by a factor of two, with the asymmetric model having 
the larger uncertainty. The wider variation in <j a in the asymmetric Gaussian model allows 
the possibility of a significantly broader distribution to smaller separations. This would 
result in a large p o pulat ion of, currently, undetected spectroscopic binaries, as predicted in 



Jeffries fc Maxtedl fl2005f ) 



The tails of the N and log(a ) distributions are correlated. Figure M displays the two- 
dimensional, marginalized posterior distribution of those two parameters for the symmetric 
Gaussian model. The tail to small values of log(ao) corresponds to the tail of the companion 
rate to larger values. This is because no current imaging survey is sensitive to companions 
at separations less than 1-2 AU (see Figure H]). The Bayesian algorithm allowed parameter 
values such that there are a large number of companions at very small separations, where no 
observations exist to constrain the overall binary fraction. This reconstructs the observed 
data, but with lower likelihood (you must move 2-3cx from the best fit values in Figure [9]) 
than the best fit model parameter set. This tail does permit t he existence of significan t 



population of very closely separated companions, as predicted by I Jeffries fc Maxtedl (120051 ). 
Further studies must be carried out to determine the frequency at separations less than a 
few tenths of an AU. While the best method for this is a spectroscopic survey for either 
radial velocity variations or double-lined systems, there are other techniques that may be 
used. These methods are widely varying and include sources wi th known paralla xes which 



are over luminous, astrometric wobbles, spectral inconsistencies ( jCruz et al.ll2004l ). and over 
luminous young cluster members. 

A further intriguing difference between the results of the two different models can be 
seen in panel c of Figure [HJ There is an auxiliary peak at o a = 0.2, which suggests there 
is a secondary solution. This solution is quite different than the two main solutions, see 
best fit parameters values in Table 1. The uncertainties, while lower for the separation 
parameters, are considerably higher for the mass ratio distribution parameters than the two 
main solutions. This solution is for a very narrow separation distribution, but a considerably 
higher binary frequency (28%). The higher binary rate is necessary to produce the correct 
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number of companions in the wings of the separation distribution to compensate for the 
narrower value of a a . Finally, the relative likelihood of the secondary solution compared 
to the primary solution is significantly different. These likelihoods are the result of the 
calculation of Equation [TT] and the difference is —100. Recall, that the likelihoods are given 
by the natural logarithm of the probability. The actual difference in probability between the 
primary and secondary solutions is e -100 , so the secondary is considerably less likely. 



5.2. Mass Ratio Distribution 



One of the most significant results is seen in panel d of Figures [7] and [HI on 7, and is 
reproduced in Figure [10J As noted earlier, there was a selection bias towards ultracool dwarf 
companions to near-equal mass, tightly bound systems. Using the correction factor from 
Equation [1] and the Bayesian analysis, I can determine if this is just the result of the bias, 
or if it is a real feature of the ultracool dwarf companion distribution. 

The most likely value of 7 for the companion mass function is ~1.8 ± 0.6, to 68% 
confidence. While this parameter is not well constrained it is different from the Salpeter 
slope (—2.35) for very ma ssive stars to 7a, the index for G, K, and early M stars (—1.05 ± 
0.15) (IReid &: Gizislll997l ) to 5a, and the best fit model for field ultracool dwarfs from A05 
(— 0.3±1.0) to 3 . 5a. V alues of 7 below 0.0 are ruled out to the 99.95% confidence level. 
Burgasser et al.l (120061 ) finds a much larger value of 7, 4.2 ± 1.0. This steeper value may 
be due to the lack of M dwarf primaries in their analysis, which concentrated on L and T 
dwarfs. The large number of late-M dwarfs in this work may have softened the slope of 
the mass ratio distribution. Thus, it may be that the mass ratio distribution continues to 
steepe n even further wi t h dec reasing primary mass. Overall, the values of 7 calculated here 
and in iBurgasser et al.l (120061 ) imply that the mass distribution of brown dwarf companions 
is fundamentally different from that of isolated ultracool dwarfs, marked by the slashed 
rectangle in Figure [10J This is also a different companion distribution from higher mass G and 
K stars, but continues the trend of lower-mass primaries (M and later) having larger fractions 
of q = 1 companions. For comparison, Figure [TT1 displays the mass ratio distributions plotted 
in Figure [1] with the best fit companion mass ratio model. 



5.3. Separation Distribution 



The most significant separation distribution result of this analysis is that the dearth of 
wide companions is a real feature of the ultracool binary distribution. As seen in panels b 
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and c of Figures [7] and [H the posterior distributions for the separation distribution is well- 
constrained with a Gaussian center of ~7.2 AU^\'] • Figure [12] is the same as Figure [2] 
but with the best fit separat ion distribution display ed as well. The binaries that appear at 
separations > ~30 AU in the lBurgasser et al.l (120071 ) data are all companions associated with 
primaries not included here and thus were excluded from this analysis. However, the overall 
frequency of wide companions (>15 AU) predicted by the best fit models is ~1— 2% (Table 1). 
This is still consistent with the relative fraction of wide companions if the recently detected 
wide, young systems are included. It should be noted that wide ultracool companions to G, 
K, and early-M stars is a different matter. 



5.4. Missing Companions? 

For the two primary best fit solutions, I determined the binary fraction missed by 
the surveys used in this work as a function of separation. Figure [13] displays the results. 
Note that the fit for the asymmetric Gaussian produces a similar number of missing close 
companions as the symmetric Gaussian. This is because there is sufficient sensitivity in the 
current surveys to see a downturn in the binary frequency to smaller separations. However, 
this is a lower limit to the missing fraction. The surveys considered here do not have any 
sensitivity to very close companions and these results may simply be caused by the lack 
observations. 

The vast majority of the missed binary fraction is for separations less than 10 AU. 
These missing companions are low-q systems that the surveys were not sensitive to and tight 
spectroscopic systems with a range of q's. The overall fraction missing is 6% — 7%. This 
means that an additional ~6% of ultracool primaries have undetected companions. This 
is perfectly consistent with observational results. The observed binary fraction is ~15%, 
so when the missing fraction is added to the observed we recovered the best fit values 
for the ultracool dwarf companion frequency. Thus, according to my predictions, current 
observational efforts have recovered 60% to 75% of ultracool dwarf binaries. 

I have also tested the possibility that there is a large reservoir of low q binary systems 
to which the surveys u sed in this work were no t sens itive. In the ana l ysis o f main sequence 



spectroscopic binaries, lHalbwachs et al.l (120031 ) and iGoldberg et al.l (12003] ) show that the 



mass ratio distribution is bi-modal. There is a peak in the distribution near q of 0.8 to 1.0, 
similar to that found here, as well as a broad peak around a q of 0.2. 



To simulate this sort of distribution, I determined the distribution of delta magnitudes 
that corresponds to a broad peak around q of 0.2. This is generally reconstructed, for the 
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.ff-band, as a Gaussian distribution centered at AMk of 6 with a width of 4 magnitudes. 
To test the detectability of this low q distribution, I created a series of AM K distributions. 
I started with the best fit symmetric Gaussian model and simply added the new low q 
extension. I normalized this extension to different additional binary frequencies, N a , in the 
range of 5% - 45%. I chose this range to correspond to, at a minimum, the 1 a upper limit 
on the value of N and, at the maximum, t o giv e near parity between the high q and low q 
peaks, as is seen in lHalbwachs et all (120031 ) and iGoldberg et all (120031 ). The new parameter 
N a is just an additional binary frequency, thus a model with N a = 10% has a total binary 
frequency of 30% (recall the baseline rate of the best fit symmetric Gaussian model of 20%). 
Figure HU displays an example of this modified delta K magnitude distribution. Here we see 
the sharp peak at small delta magnitudes from the best fit symmetric Gaussian and the new 
broad distribution of low q companions. 

I determined the number of compa nions that would have been detected at AMk > 2 



in th e surveys conducted in th e K-band (IKoerner et al.lll999l ; IClose et al.ll2003l ; ISiegler et al. 



20031 . 120051 ; lAllen et all 120071 ) by simply multiplying the window function by the revised 
model. These numbers were then compared both to the number actually detected (2) and 
the number predicted by the best fit model (1.2). The number derived for each new model 
delta magnitude model is displayed in Figure [151 At the 1 a level only the lowest value of N a 
is barely consistent with the observations, which is also consistent with the uncertainties on 
the fit of N. Thus, it is not likely that there is a significant population of low-q companions 
that have not been detected. 



5.5. Reliability of the Results 

A goodness of fit test is used on the best fit models to see if the companion model form, 
Equations [2], [3], & [51 accurately represents the data. A standard Monte Carlo technique is 
used. Five hundred realizations of each best fit companion distribution model are randomly 
drawn. These realizations are constrained to have the same number of detected companions 
as the real data sets. The likelihood of the best fit to each fake data set is calculated by 
running them through the Bayesian algorithm (Equation [TTT) . If the best fit model truly 
reflects the real data then the likelihood of the best fit of the observed data should fall near 
the peak of the fake maximum likelihoods. Figure [H] displays the results, and the best 
fit model does lie well within the likelihoods calculated from the fake data sets. However, 
both the symmetric and asymmetric Gaussian models produces similarly good fits, so we 
cannot, at this time, differentiate between the two. Thus, both models can reconstruct the 
observed data, the main difference being that the asymmetric Gaussian model may provide 
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for a substantial population of undiscovered close binary systems. 

I also conducted a similar exercise for the secondary asymmetric Gaussian solution. As 
can be seen in Figure [T71 this fit is not good. The peak of the fake likelihoods lies to the 
right of the best fit values, so that model does not provide a good representation of the data. 
Overall, the ultracool dwarf companion models do a reasonable job of reconstructing the 
data and both must be considered valid at this time. 



6. Discussion and Conclusions 



While the constraints placed by my Bayesian analysis are not tight they do reinforce 
the significance of three trends already noted in the ultracool dwarf binary distribution. 
First, the peak at q ~ 1 is a real feature. The sensitivity of the high-resolution and wide- 
field imaging surveys used in this analysis, while falling to smaller mass ratios, does not fall 
off significantly until q's of about 0.7-0.5. Therefore the statistical significance of the peak 
towards high mass ratio systems is genuine as seen in the posterior distribution for 7 (Figure 
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Second, the severe lack of companions at separations greater than ~15 — 20 AU. At 
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binary fractions of ~1%. This is a very robust result. Additionally, I tested the idea 
of an undetected spike in the mass ratio distribution at low (g~0.2), and it was shown 
that a significant population was not consistent with the observations (Figures [T41 and [T5"j) . 
Between all of the imaging surveys considered in this paper there were 361 observations for 
companions wit h separations between 20 AU and ~1000 AU and zero detections. However, 
recent work by iBilleres et al.l (120051 ) using DENIS objects has uncovered one very wide 
candidate companions in a survey of 250 potential primaries. Therefore, the rough binary 
fraction of wide companions is ~1%, which is in good agreement with my predicted results 
from the symmetric Gaussian. 



T here have been four detections of wide companions to youn g ultracool dwarfs (ILuhman 



2004al . I2005L 120061 : IChauvin et all 12004 l2005t IClose et all 120071 ). The appropriate selection 
effects for these surveys are not well quantified making it difficult if not impossible to de- 
termine a binary frequency, thus precluding the addition of these results to my Bayesian 
analysis. Nevertheless, the greater prevalence of wid e and low mass r atio s ystems among 
young objects is intriguing. As has been discussed in iBurgasser et al.l (120071 ). this is based 
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on very small number statistics (8 binaries) and, therefore, requires a sign i ficant amount of 
further study to verify. However, the overall low rate from iBilleres et all (120051 ) of ~1% is 
perfectly consistent with the outcome of my Bayesian analysis and the Gaussian separation 
distribution models which predict a similar frequency of wide companions. 

Third, the decreasing of the binary fraction with decreasing spectral type primaries. I 
predict an overall ultracool dwarf multiplicity of ~20%. This continues the decreasing trend 
from G stars, ~60%, to M stars, ~30% — 40%, to my result. My results are also consistent 
with the mos t recent em p irical estimates in the literatu re, including the ~24% overall binary 
frequency of lReid et all (bood ) and the 20% - 25% of iBasri fc Reinersl (bood ). Clearly, on 
overall observational picture is taking shape, which, as my Bayesian analysis demonstrated, 
has great statistical significance. 

There are implications for the formation mechanism models based on the conclusions of 
this work. The two empirical facts that must be accounted for are the steep mass ratio distri- 
bution that peaks at q — 1 and the ~1% frequency of wide systems. Unfortunately, few of the 
numerical simulations produce concrete predictions for the distribution of b inary systems. 



Some notable exceptions are the hierarchical triple system decay models by lUmbreit et al. 



( 120051 ). which produces a separation distribution with a peak at a projected separation of 
3 AU, as well as sharp declines to both tighter and wider systems. This qualitatively agrees 
with the best fit solutions from my analysis. The peak of the distribution is roughly half of 
what I find (3 AU instead of 7 AU). However, their distribution reproduces the few tighter 
and wider binary systems of my best fit models. Thus, the Umbriet et al. models are in 



agreement with my analysis of the data. The work of ISterzik fc Durisenl (120031 ) on dynami- 
cal interactions of small-N clusters was able to qualitatively reproduce the decreasing binary 
fraction with decreasing primary mass. This is another result that I confirm in this work. 

One model that has received a grea t deal of attention i n the last several years is that 
of proto-stellar ejection, as pioneered in iReipurth &: Clarkd (120011 ). The Bate et al. group 
( IBate. Bonnell. &: Brommll2002l . 120031 ; iBate &: Bonnel]||2005l ) has produced some of the most 
complete simulations of turbulent fragmentation and ejection of brown dwarf embryos. The 
binary fractions that they produce are 5%-8% for the separation ranges that they probe 
(> 2.5 AU). Their combined work generated ~20 binaries, only one of which was a wide 
system (20 AU). To compare these predictions to the results of this work I determined the 
binary fraction of the best fit models at separations greater than 2.5 AU. This matches the 
Bate et al. separation range, and both the symmetric and asymmetric Gaussian models yield 
binary fractions of ~14%. This is a factor of 2-3 higher than predicted by the ejection model 
simulations. Thus, the ejection hypothesis of ultracool dwarf formation must, as illustrated 
by the simulation results compared to those here, be held in doubt as it is unable to reproduce 
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the ultracool binary distribution. 



There are three major observational efforts that need to be undertaken before the ob- 
servations can be considered complete. First and foremost is the current lack of large spec- 
troscopic su rveys for companions at l ess than ~1 AU . This effort is just beginning through 
the work of Jeffries fc Maxted ( 2005 ); Joergens ( 2006 ); Basri fc Reinersj (2006). The results 



from these prelim inary surveys vary signi fi cantly, from 32 — 45% in I Jeffries fc Maxtedl (120051 ) 



to 26% ± 10% in lBasri fc Reiner j (bood l. Ijeffries fe Maxtedl fcooj) found their binary fre 



quency thro ugh statistical in f erenc e of 2-3 radial velocity measurements of many objects, 
whereas the iBasri fc Reinersj (|20061) found theirs through consistent monitoring of a small 
nearby sample. iJoergend (120061 ) examined several young objects in the Chameleon 1 star- 
forming cloud and did not find a single confir med companion. The results of this work are 
in agreement with the IBasri fc Reinerd (120061 ) findings, and demonstrate that the inferred 
Jeffries fe Maxtedl (bo05h binary fractions are likely too high by a factor of 2. Table 1 shows 
the predicted binary fractions for tight (<1 AU) systems for both the symmetric and asym- 
metric Gaussian models. I find that only an additional 3%- 1 -. 
have such tight binary systems. The disc repancy between the 



of ultracool dwarfs should 



Jeffries fc Maxtedl (120051 ) and 



the concurring results of this work and IBasri fc Reinerd (120061 ) highlights the need for a 
well-defined volume-complete sample of spectroscopic binary observations. Such a sample is 
crucial to constraining this sparsely sampled region of parameter space. 

Second, is a region of parameter space that is nearly unconstrained, wide very low-mass 
ratio companions. It is unlikely that there is a significant reservoir of unseen extremely low- 
mass companions to ultracool dwarfs (Figure [15]), but it must be explored for completeness. 
Also note that the majority of the wide companions found thus far have low mass ratios. 

Finally, high-resolution imaging surveys need to be improved to reliably reach into 
separations on the order of a few tenths of an AU. This is because the sensitivity of the 
radial velocity searches falls off at this point and the sensitivity of the current imaging 
surveys begins to fall off at around 3-4 AU. This leaves an insensitive hole in efforts to observe 
binary companions to low mass stars and brown dwarfs near where we are currently seeing 
a peak in the distribution! This must be resolved to accurately measure the distribution of 
companions. 
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Fig. 1. — Comparison of mass ratio distributio ns for G stars (jDuquennoy fc Mayor! (119911 ). 
dashed), early- M stars ([Fischer fc Marcvl (119921 ). dot-dashed), and currently known ultracool 
systems ( jBurgasser et al.l (120071 ). solid). Note that the mass ratio distribution appears to 
evolve from higher to lower mass primaries. 
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Fig. 2. — Histo grams of the separation distribution of all known companions to late-M and 
L field dwarfs (jBurgasser et al.l (120071 ) . solid), and for comparison companions to G-dwarfs 
(jDuquennoy fc Mayorl (119911 ). dashed). Note that the separation distribution of early-M 
dwarfs is very similar to that of G dwarfs, according to iFischer fc Marcyl (119921 ) and is not 
plotted separately here. This plot demonstrates the lack of wide companions as compared to 
the stellar companion population and the Gaussian shape of the ultracool dwarf separation 
distribution. 
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Fig. 3. — Left: Histogram of the distribution in spectral type of the individual target pri- 
maries of the surve ys described in Section El Right: Distanc e distribution of the surveys: 
Ground-Based AO (laoseel^lJbO-Oal: Isieerler et alJEoO^. boOfih fred asterisks ), WFPC2 HST 



jReid et al.l boOll : Idzis et al.l [2OQ3I : iBouv et alil2003l: iBurgasser et al.l l2003h (green circles) 



Keck (Koe rner et al, in prep) (blue x's), IRTF (jAllen et al.l 120071 ) (cyan squares), and HST 
NICMOS JReid et al.l bood : burgasser et al.lbood ) f purple triangles) . 



Fig. 4. — The observed window functions for the combined ultracool dwarf companion sur- 



veys in delta magnitude and projected separation in AUs. The top panel r epresents t 



re ob- 



servations of the opti cal HST surveys (jBouy et al.ll2003l : iGizis et al.ll2003l : iReid et al. 



2001 



Burgasser et al.1120031). the m iddle panel the near infrared NICMOS observations (IReid et al. 



2006 



Burgasser et a 



K-band (IClose et al. 



2006 ), and the bottom panel all the near-infrared observations in th e 



2003 



Siegler et all 1200.1 1200.4 lAUen et all 120071 : iKoerner et al.lll999l ). 



The white (black in the middle panel) crosses in both panels represent the known binaries 
from the respective surveys. The scale bars give the number of observations that corresponds 
to a particular color. 



-29- 




Fig. 5. — Comparison of the raw separation distribution model in log(a) (solid histogram) to 
the transformed projected separation (log(s)) (dotted histogram) for the symmetric Gaussian 
model (left) and the asymmetric Gaussian model (right). 
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Fig. 6. — Comparison of the raw mass ratio distribution model (left) and the transformed 
model in delta magnitude (right) in this case Am^. Also compared are two different values 
of the mass ratio distribution parameter 7, the power-law index, 7 = — 1 (solid lines) and 
7 = +1 (dashed lines). 
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Fig. 7. — Marginalized posterior distributions on the four parameters of the symmetric 
Gaussian binary companion distribution model: a) The companion rate per star at near 
equal masses, N; b) The center of the Gaussian separation distribution, log(a ); c) The 
width of the Gaussian separation distribution, er a ; d) The power-law index of the mass ratio 
distribution of companions, 7. Note the long tail on N to high values and the tail on log(ao) 
to smaller values. These tails are correlated, as seen in Figure M 
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Fig. 9. — Marginalized 2D posterior distribution on the companion fraction, N, and the 
Gaussian center parameter, log(ao) for the symmetric Gaussian model. The contours denote 
10% confidence intervals. Note how these two parameters are correlated. A smaller value of 
log(ao) requires a larger value of iV in order to recreate the correct number of detections at 
larger separations. 
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Fig. 10. — Marginalized posterior probability distribution for the power-law index 7, for the 
companion mass ratio. This is reproduced from Figure Ed. The correspon ding; range of mass 
function power-law indices derived for field objects in I Allen et al.l (120051 ) is displayed as a 
vertical slashed rectangle. 
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Fig. 11. — Similar to Figured! but also displays the Bayesian best fit mass ratio model (thick 
solid lines) for the symmetric Gaussian (left) and the asymmetric Gaussian (right). Note 
that the best fit is generally consistent with the prevalence of near-equal mass companions 
to ultracool dwarfs. 
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Fig. 12. — Similar to Figure [2j but also displays the Bayesian best fit symmetric Gaussian 
separation distribution model (solid histogram). Here we see that there is likely a modest 
population of smaller separation binaries that will be discovered with current and future 
spectroscopic surveys. 
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Fig. 13. — Predicted missing binary fraction for the symmetric Gaussian (solid histogram) 
and the asymmetric Gaussian (dotted histogram) models. The binary fraction displayed 
here is the overall, not the relative. For example, these models predict that an additional 
2.2% of ultracool primaries have undetected companions at separations of 1 AU. 
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Fig. 14. — Displays the sum of the symmetric Gaussian best fit AMk distribution plus the 
added low q test distribution. The test distribution used here is for an added companion 
fraction of 45%. The solid histogram is the sum of the two distributions while the dotted 
histograms are for the two separately. 
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Fig. 15. — Displays the number of companions that would have been detected at AM^'s 
greater than 2 (open triangles) as a function of the added binary frequency of the low q test 
distributions. Also shown are the number of companions predicted by the best fit symmetric 
Gaussian model (dashed line), the number of companions actually detected (solid line), the 
la uncertainties of that measurement (dotted rectangle). 
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Fig. 16. — Histogram of the relative Bayesian maximum likelihoods calculated using 500 
realizations of the best fit companion model for the symmetric Gaussian (left) and the 
asymmetric Gaussian (right). The x-axis represents the maximum likelihood of each model 
realization compared to the best fit maximum likelihood (zero on this scale), and the y-axis 
gives the number of realizations at each relative likelihood. The maximum likelihood for 
each realizations was compared to that derived for the best fit model. If the model is a true 
representation of the data then these histograms should peak near zero, which they do. 
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Fig. 17. — Similar to Figure fT6l but for the secondary solution to the asymmetric Gaussian 
model. This demonstrates that this secondary solution is not a good fit to the data. The 
likelihood of that solution (at zero on this scale) sits well off of the peak. 
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Table 1. 



Parameter Symmetric Asymmetric Secondary 

Gaussian Gaussian Solution 



log(a ) [log(AU)} 


0.861^ 


o.86ig;Si 


0.86 ±0.06 


a a [hg(AU)} 


0.28 ±0.04 


0.30 ±0.10 


0.20 ±0.02 


N 


20% ± 4% 


22%1J* 


28% ±11% 


7 


1 8+ - 4 
L °-0.6 


1.6 ±0.6 


1.4 ± 1.4 


% Comp. > 15 AU 


1.1% 


1.3% 


NA 


% Tight Binaries 


2.9% 


3.5% 


NA 



